Astronomy & Astrophysics manuscript no. HS201 1-ArXiv 


©ESO 2011 


October 25,2011 





Star formation near an obscured AGN 

Variations in tlie initial mass function 

S. Hocuk^ and M. Spaans' 



Kapteyn Astronomical Institute, University of Groningen, P. O. Box 800, 9700 AV Groningen 
e-mail: seyit@astro . rug . nl , spaansOastro . rug . nl 



o 
o 

o 

(N 

u 

o 



> 
cn 

o 



X 



Received October 25, 2011 



ABSTRACT 



The conditions that affect the formation of stars in radiatively and mechanically active environments are quite different than the 
conditions that apply to our local interstellar neighborhood. In such galactic environments, a variety of feedback processes can play 
a significant role in shaping the initial mass function (IMF). Here, we present a numerical study on the effects of an accreting black 
hole and the influence of nearby massive stars to a collapsing, 800 Mq, molecular cloud at 10 pc distance from the black hole. Our 
work focusses on the star-forming ISM in the centers of (U)LIRGS. We therefore assume that this region is enshrouded by gas and 
dust and that most of the UV and soft X-ray radiation from the BLR is attenuated along the line of sight to the model cloud. We then 
parametrize and study radiative feedback effects of hard X-rays emanating from the black hole broad line region, increased cosmic ray 
rates due to supemovae in starbursts, and strong UV radiation produced by nearby massive stars. We also investigate the importance 
of shear from the supermassive, 10^-10** Mq, black hole as the star-forming cloud orbits around it. A grid of 42 models is created and 
calculated with the hydrodynamical code FLASH. We find that thermal pressure from X-rays compresses the cloud, which induces a 
high star-formation rate early on, but reduces the overall star-formation efficiency to about 7% due to gas depletion by evaporation. 
We see that the turn-over mass of the IMF increases up to a factor of 2.3, M[„r„ = 1-1.5 Mq, for the model with the highest X-ray 
flux (160 erg s"' cm"-), while the high-mass slope of the IMF becomes F > -1 (Fsaipctcr = -1-35). This results in more high-mass 
stars and a non-Salpeter initial mass function. Cosmic rays penetrate deeply into the cloud and increase the gas temperature to about 
50 K for rates that are roughly 100 times Galactic and 200 K for 3000 times Galactic, which leads to a reduced formation efficiency 
of low-mass stars. While the shape of the mass function is preserved, high cosmic ray rates increase the average mass of stars, thereby 
shifting the turn-over mass to higher values, i.e., up to several solar masses. Due to this process, the onset of star formation is also 
delayed. We find that UV radiation plays only a minor role. Since UV photons cannot penetrate a dense, n > 10^ cm"^, cloud deep 
enough, they only affect the late time accretion by heating the medium where the cloud is embedded in. When we increase the black 
hole mass, for a cloud that is at 10 pc distance, the turbulence caused by shearing effects reduces the star-formation efficiency slightly. 
Furthermore, shear weakens the effect of the other parameters on the slope of the IMF as well as the turn-over mass. The run with 
the most massive black hole, however, causes so much shear that the hydrodynamics is completely dominated by this effect and it 
severely inhibits star formation. We conclude that the initial mass function inside active galaxies is different than the one obtained 
from local environments. We also find that the combined effects of X-rays, cosmic rays, UV, and shear tend to drive toward a less 
pronounced deviation from a Salpeter IMF. 

Key words. Stars: formation - Stars: mass function - X-rays: ISM - Cosmic rays - Radiative transfer - Methods: numerical 



1. Introduction 

Star formation in extreme environments can be quite different 
than the formation of most stars in the Universe. In the inner kpc 
of galaxies, molecular clouds are exposed to in tense radiation 
from active galactic nuclei (AGN ) or starbursts dMaloney et all 
19961: iMe iierink & Spaaiisl 120051: iPerez-Beaupuits et all 120091 



van derWerf et al.l l2010h . Very close to the AGN, < 0.1 pc, gas 



collects into a massive AGN disc over som e time, and as the disc 
becomes unstable, star s are able to form dPaumard et al.ll2006t 
iNavakshin et al.l l2007h . Slightly further away from the black 
hole, 1-100 pc, conditions are somewhat less extreme as the ra- 
diation is strongly attenuated by large columns of gas and dust. 
Such environments are typical of obscured AGN, as formed in 
(ultra-)luminous infrared gal axies ((U)LIRGS), with obscuring 
coluir i ns of 10^^-10^^-^ cm~ ^ (lAaltol2005l;lPerez-Beaupuits et all 
120071: |L oenen et al ] I2008h . These regions have a strong im- 
pact on the initial phases of cloud evolution so that the fi- 
nal mass of stars or their formation efficiencies might drasti- 
cally change. However, observing star formation in extreme en- 



vironments is difficult. Results usually rely on indirect meth- 
ods and are therefore ofte n debated. The regions near AGN 
are al so generally obscured (ISpoon et al .120071: iHocuk & Barthell 
I20T0I) . Predictions based on models and numerical simula- 
tions can aid observations to further our understanding of 
star formation. A good amount of numerical work has been 
done focussing on mechanical and radiative f eedback effects 
in active galaxies a n d star-forming r egions (iKlessen 



KlessenetaTI 12005 at fBonnell & Ric3 120081: [Alexander etal 



2008';'Wada"2008';'Hobbs & Navakshin"2009'; 'Wada et al."2009 
Krumholz et al. 2010; Bate 2010; Hocuk& Scaans 2010b^; 
Perez-BeauDuits et al. 201 1; Clar k et al.ll201 ll Latif et al.ll201 1 



Avkutalp & Spaans..2011; Alia et al.l 1201 lb . Supported by nu- 



merical simulations, it is often thought that in these active re- 
gions, due to environmental conditions and feedback effects, the 
IMF should be different than the proposed universal mass func- 
tion (Klessen et al. 200 71 INavakshin et al.l200"7tlHsu et al.l2010l: 
iKmmholz et al.,i20ia:iGirichidis et al.ll201 ll) . 



2 



S. Hocuk & M. Spaans: Star formation in AGN 



The initial mass function is of fundamental importance for 
many areas of astrophysics. It is observed to behave like a power- 
law with a high-mass end that is well defined. The IMF is de- 
scribed as 



dN 
dM 



oc M 



dlogN 
dlogM 



-a H- 1 = F, 



(1) 



with N the number of stars in a range of mass dM, a the 
power-law index, and F the slope a bove th e characteristic mass 
of -0.3-0.5 Mq. First proposed bv ISalpe ter (1955), a plethora 
of observations has led astronomers to believe that the shape 
might be universal in nature. Other astronomers have refined 
the shape of the distribut ion by especially improving o n the 
low-mass end of the IMF ( iMiller & Scalolll979t iKroupalllOOlt 
[Chabrier 2003). Until observations of extragalactic origin started 
to show variations in the IMF, most of the studies that claimed 
universality were coming from observations from our local 
neighborhood. Hints for deviation came fro m measurements of 
abundance patterns in extragalactic bulges (iBallero et al.ll2007l 
l2008h . enhancement of far infra-red luminosities in interacting 



galaxy systems dBrassington et al.| 2007h. mass-to-light ratio s 
of ultra-compact dwarf galaxies dPabringhausen et al. l2009h. 
and many others ( Baugh et"al] 20051: Parra et alT 2007 : iDayel 



,2008VWilkins et al 2008; Van Dokkum 2008;'Elmegreen 2009*). 
vari DoIJmm &; Conroy (2PJ,0j 2D11) found further evidence for 
variations in the IMF at the low-mass end from Nal and FeH 
band spectra in luminous elliptical galaxies. Still, strong evi- 
dence has yet to emerge. 

Although the initial mass function may seem, theoretically, 
to be univers al over a relatively wide range of enyironmen- 
tal conditions (IClark et al.ll200 9i iBastian et al 1 1201 0HKmmhoi3 



1201 1*). perha ps even in sensitive to small changes in metalUc- 
ity (M yers et aP |20I ih and the Jeans mass (Elmegre en et alJ 
120081) . there are conditions that are far more extreme than 
the ones disc ussed in these papers like in radiation domi- 
nated regions (Meiierink et al.ll2005t Poelman & SDaansll2006t 



Ivan der Werf et alJl20l tf. 'Meiierink e t al. |20I lb and co smic ray 
dominated regions (CRDRs, Papadop oulos et a0201 lb . Besides 
the importance of the thermodynamics for the Jeans mass, Mj, 
and thus the IMF, where Mj is proportional to p"'^-T^^^, the 
change in the equation of state is also essential. Assuming ideal 
gas conditions with a polytropic equation of state, P oc p^, where 
7=1-1- dlog(T)/dlog(p), the softness of y plays a major role, at 
a very early stage, in the f ragmentation properties and the mass 
scale of unstable clouds dSpaans & Siljd 120001: iLi et alJ l2003t 
iKIessen et al ] l2005bHJappsen et alJl2005h . 

The nuclei of active galaxies like Arp 220, Markarian 
231, and e ven our Gala cti c center show signs of unusual star 
forma t ion (iFiged l2005allbl: iPaumard et alJ l2006t iKlessen et al] 



2007|: lEspinoza etal.l 120091: 'Elmegree n! 120091: 



201C ;|Mi.tsushita et al. 2009; Meiierink et alj20l U 



I et alJ 
I et alj 



Bartko i 
Martin et al.i 

201 II) . In fact, even in the inner parsec of our Galaxy, i.e., Sgr 
A* and in M31, youn g stars are found a t distances on the or- 
der of - 0.0 3 - 0.3 pc (iGenzel et alJl2003l:lPaumard et al.ll2006t 
lLevinll2007l) . All of the aforementioned places harbor a mas- 
sive black hole. One can imagine that the conditions close to 
the black hole can indeed become quite extreme. Aside from 
strong gravity, accreting material onto a black hole will produce 
strong X-ray radiation (1-100 KeV). The dynamics of molec- 
ular clouds will be significantly affected by the irradiation of 
X-rays in X-ray dorn inated regions (XDRs, J^epp & Dalgarnoi 
Ima iMalonev et al.1 [19961) . On the other hand, in starbursts, 
where star-formation rates can be a few hundred to a thousand 



solar masses per year dSmail et al.lll997t iHughes et al.lll998h . 
UV radiation (6-13.6 eV) from O and B stars can be a signifi- 
cant presence and dominate the ra diation field in these so-calle d 
photon dominated regions (PDRs, iHoUenbach & Tielensl[T999l) . 
However, where the gas is shielded from UV radiation, cosmic 
rays, c reated in supernova remnants or from winds in OB associ- 
ations dBinns et al.ll2008h . will dictate the (minimum) tempera- 
ture of the syste m, with energy densit i es of up to a few thousand 
times our galaxy (Patjadopoulosl20 lot iPapadopoulos et alj20I II : 



iMeiierink et al.l2bl ih . All of these environments have extremely 
different star-forming conditions, but will the stars that form out 
of them be much different? 

In radiation dominated regions, the chemi stry and thermal 
balance are determined by the radiation field dSchleicher et al.l 
I2OIOI) . X-rays, in active galactic nuclei, are usually the dominant 
source for the excitation and chemistry of the inner disk out to 
a radius of ~160 pc, while UV radiation dominates the thermal 
balance in extreme starbursts and is genera lly important a bit 
further away from the accreting black hole dvan der Werf et all 
I2OIOI) . But there is compelling evidence that there is a strong 
link between AGN and starbursts (Sc oville 2004.). The question 
remains, how strongly star formation is affected by these extreme 
environments. 

In an earlier numerical study, we showed that when the X-ray 
flux is as high as 160 erg s"' cm"^, the stellar initial mass func- 
tion of an 800 solar mas s molecular cloud becomes top-heavy 
dHocuk&Spaansll2010ah . This was a case where a molecular 
cloud orbiting at 10 pc from a 10^ M© black hole under the im- 
pact of X-rays was compared to a cloud with the same conditions 
but in an X-ray free environment. This numerical study showed 
that under extreme conditions, the evolution of a molecular cloud 
and its stellar mass function will change, but it did not give in- 
sight into the quantitative details. Here, we present a parameter 
study on the influence of external radiation (X-rays, cosmic rays, 
and UV) and black hole shear on the IMF and the star-formation 
efficiency (SEE). In section |2] we introduce the numerical code 
FLASH and describe our additions to it. In section [3] we define 
the cloud models for all ambient conditions considered in the 
parameter study. We then present in section |4] our results on the 
effects of each condition for the evolution of the model clouds 
and show their phase diagrams, star-formation efficiencies, and 
initial mass functions. Finally, in section |5] we discuss the dif- 
ferences and similarities of these results in detail and present our 
conclusions. 



2. Computational method 

2. 1 . The numerical code 

The calculations in this wo rk have been done using the hydro- 
dynamical code FLASH 3 (iDubev et al.ll2009l) . For this study, 
we use the directionally split piec ewise-parabolic metho d (PPM ) 
which is described in detail in IColella & Woodward! dl984 . 
FLASH is well suited to handle these kinds of calculations as 
it is an adaptive mesh code and one that handles contact discon- 
tinuities very well. FLASH is provided with many and exten- 
sively tested modules that encompass a broad range of physics. 
Our simulation code is equipped with thermodynamics, hydro- 
dynamics, (self-)gravity, multi-species, particles, and shocks, all 
from the standard modules of FLASH, as well as sink particles, 
radiative transfer, multi-scale turbulence, and refinement crite- 
ria (based on Jeans length and particles) that were added by us. 
The non-standard additio ns are explained in furthe r detail in the 
following sections and in lHocuk&Spaansld2010ah . 
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2.2. Refinement criteria 

When one does a parameter study and has to perform many nu- 
merical simulations, saving time becomes crucial. In order to 
achieve this without suffering loss of quality, we made use of the 
adaptiveness of the FLASH code and wrote two independent re- 
finement criteria that served our purpose. The simpler one of the 
two is based on sink particles. Since every sink particle accretes 
matter, and this matter can only be followed to within the sink 
particle's accretion radius, it is best to have the highest possible 
resolution here in order to resolve the aff'ected volume properly. 
Particles can only be created in regions that have the highest grid 
resolution, however, they can move through the grid to unrefined 
regions. To this end, we simply say that wherever a sink particle 
is located, the grid must be refined to its maximum. 

The second criterion is based on the Jeans length. In or- 
der to avoid numer ical effects, like artificial fragmentation 
dTruelove et al.|[T997l) . it is necessary to resolve the Jeans length, 
Aj, in the simulation by at least 4 cells. As > 4ax. Here, Ax 
is the size of the grid cell which is dependent on the resolu- 
tion. In our experience, however, it is also likely that fragmen- 
tation can be artificially inhibited at even higher resolutions, up 
to twice the Truelove criterion, that is, 8 cells. In the presence 
of magnetic fields, even h igher resolution constraints are found 
by iFederrath et al.l (1201 ih . However, we do not consider mag- 
netic field effects in this study. Therefore, we have chosen to 
resolve the Jeans length in our simulations by at least 10 grid 
cells, Aes.gas = 10. This refinement criterion can be rewritten in 
the form of a density threshold and refines the grid if the follow- 
ing condition is met: 



M, 



|7r(^resAx)3 ' 



(2) 



with the Jeans mass, Mj, defined as: 



i 



Mj = -7r^> = p-i - p Ti n-5 (M^), 



(3) 



as taken from lFrieswiik et al] ( l2007h . Using this, Eq. ^ can be 
further reduced to: 



p> 1.455 X 10' 



jU(A-esAx)2' 



(4) 



where fi is the mean molecular weight, G is the gravitational 
constant, and T is the gas temperature. 

Since we have a body in orbit, the Jeans refinement criterion 
will automatically follow the cloud in motion and increase the 
resolution of grid cells whenever required. The border resolution 
between the refinement levels is set-up in such a way that there 
is no sharp transition between the minimum and the maximum 
refinement. In order to be efficient, we also de-refine the regions 
where the Jeans length has stretched beyond 25 times the grid 
resolution. One can imagine this as a region where the cloud has 
just passed through. 

2.3. Multi-scale turbulence 

Turbulence is an important aspect of star formation. Never is the 
interstellar matter from which stars form fully homogeneous or 
kinematically quiescent. Typical velocity dispersions are found 
to be o n the order of 1 km/s in most of the regions in our 
galaxy (lLarsonlll98ll: iFalgarone et"ani200ll: ICaselli et al.ll200"2h 



and scale according to a power-law. It should be no surprise that 
gaseous clouds within active regions are more turbulent. Typical 
FWHM values are of the order of 5 km/s IPerez-Beaupui ts et al.l 
120091: lOstriker & Shettvll20Tlh . 

Numerical simulations have shown that turbulent strength, 
scaling, or type (compressible or solenoidal) can indeed be quite 
important an d affect the re sults significantly ( IFederrath et al.l 
2009; Giric hidis et al.ll20in) . It is thus important to incorporate 
turbulence into the numerical code in a proper way. In all our 
simulations, we implement the turbulence using a Larson power- 
law, with a power spectrum P( k) oc k''^ and thus Av oc where 
k PC is the sc ale length (iLarsonl Il98lt iMvers & Gammid 
ll999l : lHever & Brunt 2004) . This is the predicted and observed 
behavior for compressible fluids. Since the grids in hydrodynam- 
ical simulations are discretized, and in block-structured grids the 
cell sizes usually increase with a factor of two for each resolution 
increment, the super-posed velocity decreases with the square- 
root of two for each higher level. In our case, the largest scale 
that we apply the multi-scale turbulence on is that of the cloud 
radius. Otherwise, the cloud as a whole would obtain a random 
motion. The smallest scale on which the turbulence is injected 
is determined by the maximum resolution at runtime, which is 
Ax = 1.76 X 10'^ cm in our runs. 



2.4. The radiative transfer method 

Radiative transfer for PDRs and XDRs is handled through a ra- 
diation dominated region code wri tten by Meiierink & SpaansI 
(20()5l). with additional de tails in iMeiierink et alj (120071) and 
ISpaans & Meiierinfl (l2008h . Pre-computed tables for gas tem- 
perature and chemical abundances are obtained from this 
code given an input of number density [cm^^], radiation flux 



[erg s" cm ] for X-rays and UV, column density [cm" ], and 
metallicity. We choose solar abundances for all simulations. 
For the regions dominated by radiation, the code calculates all 
the heating processes (photo-ionization, yielding non-thermal 
electrons, FUV pumping followed by collisional de-excitation), 
cooling processes (atomic fine-structure and semi-forbidden 
lines), and molecular transitions (CO, H2, H2O, OH, and CH). 
Cosmic rays, dust-gas coupling, and secondary effects from X- 
rays like internal UV, are treated as well. These tables are in- 
corporated into FLASH. For the cloud models in the absence of 
radiation, we use isothermal conditions with an equation of state 
of the form P oc p (Ho cuk & Spaan s 2010aii'b|). 

We use a ray-tracing method to find the column densities dur- 
ing the simulation. At each time-step, the algorithm searches the 
grid and sums up the densities of each cell lying along the line of 
sight from the radiation source in order to find the total column. 
The cells are selected by determining whether the two angles of 
the cell edges with respect to the radiation source accommodate 
the angle of the target cell. Each cell is weighted with the length 
of the ray that passes through it. Once the column density is 
found, together with the density and the flux, the corresponding 
temperature is taken from t he tables and used to update t he vari- 
ables in the simulation. See lPerez-Beaupuits et alJ(l201 ih for ad- 
ditional details. In order to gain speed, the algorithm makes use 
of the block-structured mesh of FLASH by first searching the 
blocks that lie within the line of sight of the source, thereby only 
looking into the cells of those blocks. This makes the ray-tracing 
in effect about 500 times faster This gain is welcome, since, by 
far, the most computational time is normally spent on finding 
the column densities. Additional speed is gained by not apply- 
ing the ray-tracing algorithm to the whole grid at every timestep. 
We prioritize the regions that have densities above 100 cm"-' and 
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update them regularly, but the lower density regions are only up- 
dated occasionally. These regions are of not much importance 
since stars cannot form here. Besides, the temperatures at low 
densities do not decrease quickly, since the cooling time scales 
as tcooi <^ 1/n (optically thin and sub-thermal), and regular up- 
dates are not necessary. 

The X-ray flux is an E"" power-law between 1 and 100 keV. 
X-ray scattering is not very important, but is nonetheless treated 
in the XDR-code. Inverse Compton heating is not included, be- 
cause our focus lies on the molecular gas (<1000 K) that is 
mostly h eated by ionization of H a nd H2, and by recom- 
bination dMeiierink & Spaa7iil2005h . A uniform background of 
cosmic rays prevents the temperature from dropping below 10 
K. For this, a minimum cosmic ray ionization rat e typical for the 
Milky Way, fcR.Gai = 5 x 10"'^ s"', is assumed (ISpaans & SiUd 
I2OOOI) . The UV flux enjoys en ergies betwee n 6 and 13.6 eV and 
follows the Habing spectrum ( lHabing|ll968h . 

2.5. Sink particles 

Sink particles are necessary if one wants to obtain a high 
dynamic range in density. These particles represent compact 
(proto-stellar) objects that are indivisible but can gain mass 
by accreting or merging. It is computationally very expensive 
to keep increasing the grid resolution in order to follow a 
molecular cloud collapse up to proto-stellar densities. We 
have created a sink particle algorithm for this purpose. There 
are several criteria that we check to establish irreversible 
collapse. When this stage is reached, one can stop following the 
collapse an d can m ake a transition to sin k particl es. Following 
iBate et al.l (Il995h and iFederrath etaP OOlOh . with only 
small differences, we determine the point-of-no-return of a grid 
cell within a volume that is defined by an accretion radius race as: 



1. The grid cell is about to violate the Jeans criterion, Eq|4] 

2. The grid cell and its neighbors are at the highest level of 
refinement. 

3. The grid cell has the deepest gravitational potential of all the 
cells within the volume. 

4. The divergence on each axis of the grid cell is negative, that 
is, dvi/di < 0, with i = {x,y, z), such that V • v < 0. 

5. The volume within the accretion radius is gravitationally 
bound, Egi-av + Eth + Ekin < 0. 

6. The volume within the accretion radius is Jeans unstable, 
Egi-av + 2Eth < and thus M > Mj. 

7. There are no other sink particles occupying the same cell. 

At the moment when a sink particle is created, material from 
within an accretion radius is taken away from the cells and put 
into the particle. We make sure that the total mass and momen- 
tum are conserved. The accretion radius is ideally smaller than 
the sonic radius, Cjdt, but one also wants to resolve this region 
well. On the other hand, one wants to keep the radius small, since 
one changes the physics within this radius. In practice we found 
that 2 cells was a reasonable compromise between these com- 
peting demands, also see Krumholz et al. (2004); Federrath et all 
(I2OIOI) . Our accretion radius is race = 3.5 x 10'^ cm and repro- 
duces the observed Chabrier IMF for runs MOl and M04, see 
section m 

The mass a sink particle obtains at creation is determined by 
a density threshold that is set by the maximum resolution. In this 
case, w e follow the same res olution criterion as previously men- 
tioned (iTruelove et al.l[T997h . but let the code run as long as pos- 



sible on gas dynamics before making the transition to sink par- 
ticles. Using Eqs.|2]and|4l and taking Aes.sink - 6, we determine 
how much material within race from the target cell is in excess 
of the threshold density and add this to the mass of the particle. 
Sink particles continue to accrete matter after they are created. 
Accretion onto particles is handled by a Bondi-Hoyle type of 
accretion /Bondi & Hovle"l 944l: lBondilll952l: iRuffert & Arnetj 
1994; Krumholz etjiL 2004). This kind of accretion applies to a 
homogeneous flow of matter. Accretion increases with protostel- 
lar mass M, but drops with decreasing ambient density p^a and 
higher Mach numbers Ai. The accretion rate is given by 



M = Anp^G-M^ 



(5) 



where Cm and Vco are the sound speed and the velocity of the 
gas far from the sink particle. Since there is no obvious choice 
for Coo and Vco in an inhomogeneous environment, we take the 
average value of the cells inside race as an alternative. A is a non- 
dimensional parameter that depends on the equation of state of 
the gas and is on the order of unity. Throughout this work, we 
have adopted the v alue for an isothermal gas A - exp(3/2)/4 ^ 

1.12 (lBondilll952h . 

Using the effective radius as first described bv lBondil(fT95l 



Tbh 



GM 

^00 ^ • 00 



Eq. |5]can be further reduced to 



M = 47rpcorgjjCco + M?, 
where poo is defined as 



P 

Poo - —■ 

a 



(6) 



(7) 



(8) 



Here, p is the mean density within the accretion radius and a is 
the density profile that depends on the cell size and the Bondi- 
radius. For the density profile, we use an exponential of the form 
a - exp(rBH/l-2Ax). We found this expression to behave well 
in the regime Ax > rBH, where we usually are in, and equals 
to unity when Ax » reH- The factor of 1.2 is adopted from 
iKrumholz et~an(l2004 and should give good results in the range 
where Ax ~ rsH- 

Stars are often found in binaries, however, we cannot resolve 
binaries in our numerical code as our resolution is on the order 
of 2000 AU. Instead, we allow them to merge. These mergers 
do not affect our conclusions, see section |5] Binary formation 
or mergers should occur more frequently as the system virial- 
izes. Sink particles are eligible to merge when they are within 
each other's gravitational pull. We let an algorithm check for 
three conditions between every two particles, and when they pass 
these, we allow them to merge. This happens when: 

1. The velocity difference between two particles, i and j, is less 
than the escape velocity between them, AVy < -^2GMi/Arij, 
with Mi > Mj. 

2. The merging time is shorter than the hydrodynamical 
timestep, jtA?^^/ ^J8G(M[+M0 < dthydro- 

3. The forces of the other bodies, Ftot, are no longer significant, 
i.e., Ftot < 0.05 X Fy, where Fy = GMiMj/Ar^. 
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3. Models and initial conditions 

3.1. A grid of models 

For this work, a grid that consists of 42 models is computed 
which covers a range in the 4 parameters that we investigate. 
As stated before, the parameters that we vary are the X-ray 
flux, the cosmic ray rate, black hole mass, and the UV flux. 
A 10^ M0 black hole accreting at 10% Eddington and at a 
distance of 10 pc, would radiate with a total flux of pEdd = 
10"* (Mbh/lO^Mo) erg s"' cm"^. We assume, since we are in- 
terested in (U)LIRGS, that there is a large absorbing column 
with Tikev ~ 5 and tuv S 50 between our single star-forming 
cloud and the black hole BLR. Varying the black hole mass Mbh 
or its distance from the model cloud dbh, gives the same insight 
into the effects of gravitational shear, where the largest velocity 
difference Av across the cloud depends on Mbh/dj^^ for a given 
model cloud size much smaller than dbh. We have chosen to fix 
the cloud distance to the black hole and vary the shear using the 
black hole mass. 

We study four X-ray fluxes by varying the Eddington rates, 
these are: 0, 5.1, 28, 160 erg s"' cm"^, with roughly a fac- 
tor of 5 between them. X-ray fluxes much higher than 100 erg 
s"' cm"^ would completely ionize any molecular cloud of < 
10^ cm and inhibit star formation. We take three black hole 
masses: 10^, 10^, and 10^ Mq, which, for a fixed distance of 
dbh - 10 pc, represents strong shear, medium shear, and negligi- 
ble shear, respectively. Furthermore, we investigate three differ- 
ent cosmic ray rates: 1, 100, and 300 times Galactic (^cr Gai - 
5 X 10 '^ s-\ ISpaans&Sil5l2000l) . Finally, we use two UV 
fluxes in our simulations for wh ich we take: and 10^^ Go, with 
Go = 1.6 X 10""^ erg s^' cm 2 (lHabinglll968h . For each param- 
eter that we change, we keep all other conditions fixed. Almost 
all possible combinations between these parameters, with the ex- 
ception of UV, is modeled in this study. Our fiducial model MOl 
has similar conditions to those of the Milky Way, albeit, with 
added gravitational shear In this respect, model M04, which has 
negligible shear, is more closely related to the conditions of the 
solar neighbourhood. Table [T] shows the parameter details of all 
models. 



3.2. The initial conditions 

We created a simulation box of size 24-' pc-' with outflow bound- 
aries in each direction and isolated in terms of gravity. In here, 
we put a typical molecular cloud, or clump as some prefer, with 
a uniform number density of 10^ cm""* and a size of rdoud = 0.33 
pc in spherical radius. With a mean molecular weight of fi - 

2.3, the total mass of the cloud amounts to 800 solar masses. 
The rest of the medium is filled with gas that has a uniform 
density of 10 cm"-'. The total gas mass of the simulation box 
is 8000 solar masses. At the center of the box, we put a point 
particle with mass Mbh that represents the black hole, where 
Mbh - [10'', 10'', 10^] Mq. The temperature of the gas is ini- 
tialized as 10 K, but depending on the model we expose the 
cloud, immediately after the simulation starts, to external radi- 
ation. Since these radiative processes are fast, with respect to 
the hydrodynamics of the simulation, the temperature changes 
quickly, within 10'^ s {- one timestep), after initialization. The 
radiation source is either an X-ray emitter (accreting black hole), 
a uniform background of cosmic rays (mainly supernova rem- 
nants), and/or an isotropic UV radiation field (nearby massive 
stars). 



An initially random, divergence-free turbulent velocity field 
is applied to the molecular cloud with a characteristic FWHM 
of 5 km/s that agrees well with c louds observed in active re- 
gions dPerez-Beaupuits et alJl2009l) . For the isothermal runs, the 
sound speed of the cloud is c, - 0.19 km/s (for T=10 K) and 
can go up to ~5 km/s when radiation impinges on the cloud. 
There are supersonic flows with Mach numbers of up to 25. We 
apply the turbulence over all scales with a power spectrum of 
P(k) oc k as mention ed earlier, following the empirical laws 
for compressible flui ds (lLarsonlll98lt iMvers & Gammielll999l : 
iHever & Brun3l2004l) . The turbulence in this work is not driven. 
Still, it can remain strong throughout a simulation due to gravi- 
tational instabilities or shear induced by the black hole, and does 
so for the larger Mbh runs. 

The simulations start with a cloud at 10 pc distance from the 
black hole that is in a stable Keplerian orbit. The orbital time 
is on the order of 10^ x (10^ Mg/Mbh)''" yr, which is larger 
than the cloud free-fall time in any of the models. Shear due to 
the black hole gravity, creates a maximum velocity difference of 
about AVsheai- = 2.22 X (Mbh/lO' Mq)''- km/s across the cloud. 
This follows from 



Av.he.= y^3;— , (9) 

where rdoud « dbh. The shearing time, t^hear = 2rcioud/AVsheai-, is 
about 1-10 times larger than the cloud free-fall time and gravita- 
tionally bound (roughly) spherical clouds are thus likely to exist 
at densities of ~10^ cm"-'. 

The maximum grid resolution that we allow for in any sim- 
ulation is 4096^* cells. For a box of size 24 pc, the maximum 
spatial resolution becomes Ax - 1.76 x 10'^ cm. All the sim- 
ulations are evolved up to 3 free-fall times, where one free-fall 
time, given the initial conditions, is tff - 10^ years. This is taken 
as the basic time unit throughout this work. 

4. Results 

We divide the results into four sections. In each section we in- 
dividually evaluate the effects of X-rays, cosmic rays, UV, and 
shear, while keeping every other parameter fixed. In Figs. |6] 
to [171 we first show the results of all parameter variations at 
once, as this gives a better overview. Figs. |6] to |9] display the 
temperature-density phase diagrams. Figs. [10] to [13] the star- 
formation efficiencies, and Figs. [14] to [17] the initial mass func- 
tions. 

For each phase diagram in Figs. [6] to [9] we plot the number 
density versus the temperature for one moment in time, which is 
at 1 tjf . The diagram is subsequently gridded into 75^ cells and 
the weighted masses of all the points within each cell is summed 
up and depicted in color Note that the isothermal conditions al- 
ways yield a flat profile. 

The star-formation efliciencies in Figs.[T0]to[T3]are displayed 
by plotting the ratio of the total sink particle mass over the to- 
tal initial gas mass (8000 M©), SEE = Msink/Mtotah against time. 
The efficiencies are plotted as they are, i.e., no fitting is involved. 
Generally between 1 and 2 tff, the star-formation efficiencies flat- 
ten out. This is due to the depletion of high density gas. At this 
stage, star formation is almost completely halted and accretion 
is reduced to a minimum. 

Each IMF plot in Figs. [14] to [TT] is constructed by taking 
all the sink-particles of 61 different cloud evolution times of 
one model and by normalizing them. For this, we took the time 
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Table 1. Parameter details for each model 



Model 


Fx 


Mbh 




Fuv 


Model 


Fx 


Mbh 




Fuv 




r -1 -2n 

[erg s cm J 


[Mq] 


[xGalactic] 


[Go] 




[erg s cm -J 


[Mo] 


[xGalactic] 


[Go] 


MOl 





TTTT 

10^ 


1 





M22 


160 


10* 


1 





M02 





10^ 


100 





M23 


160 


10* 


100 





M03 





10' 


3000 





M24 


160 


10* 


3000 





M04 





10'' 


1 





M25 





10" 


1 





M05 





lO"" 


100 





M26 





10" 


100 





M06 





10* 


3000 





M27 





10" 


3000 





M07 


5.1 


10^ 


1 





M28 


160 


1 rvS 

10" 


1 





M08 


5.1 


10' 


100 





M29 


160 


10 


100 





M09 


5.1 


10' 


3000 





M30 


160 


10" 


3000 





MIO 


5.1 


10* 


1 





M31 





10' 


1 


10^^ 


Mil 


5.1 


10* 


100 





M32 





10' 


100 


10- 


M12 


5.1 


10* 


3000 





M33 





10' 


3000 


10^^ 


M13 


28 


10' 


1 





M34 


5.1 


1 f\l 

10' 


1 


10- 


M14 


28 


10^ 


100 





M35 


5.1 


10' 


100 


10-^^ 


M15 


28 


10^ 


3000 





M36 


5.1 


10' 


3000 


10^^ 


M16 


28 


10* 


1 





M37 


28 


10' 


1 


10^5 


M17 


28 


10* 


100 





M38 


28 


10' 


100 


10^^ 


M18 


28 


10* 


3000 





M39 


28 


10' 


3000 


10^5 


M19 


160 


10^ 


1 





M40 


160 


10' 


1 


10^^ 


M20 


160 


10^ 


100 





M41 


160 


10' 


100 




M21 


160 


10' 


3000 





M42 


160 


10' 


3000 





frames between 1 tff and 3 tff, with a time resolution of 1/30 
t)f . The sink particle masses are then plotted as a logarithmi- 
cally binned histogram, with a fixed number of 16 bins, of stel- 
lar mass versus number We compared these IMFs with those 
corresponding to the time frames between 1 tff and 2 tff, and 
found that there was little difference between them. There are 4 
colored lines plotted in each of the sub-figures. The green line 
represents the Salpeter IMF with a slope of Fsai = -1.35 and 
is shown for comparison purposes only. The blue line displays 
the Chabrier IMF as fitted to our fiducial case, model MOl, an 
isothermal simulation at t = tff (this is a single snapshot in time) 
with Mbh - 10^ Mo,Fx - 0, Fyy ^ 0, and fcR = 1 x Galactic. 
See also iHocuk & SpaansI (l2010ah . This curve is a lognormal 
function up to 1 Mq with a power-law tail of Fcha = -1.30 
(IChabrieiil2003h . The red solid line gives the best (least squares) 
lognormal fit, whereas the purple triple-dotted line shows the 
best (least squares) power-law fit above the turn-over mass. The 
turn-over mass is not fixed and can change for each model. We 
determine the turn-over mass by finding the position with the 
lowest absolute derivative (slope) |dlogN/dlogM| — > from the 
lognormal fit. 

4. 7 . Effects of X-rays 

X-rays heat up the lower density unshielded gas and can do so 
up to 6000 K at the highest flux. Hard X-rays (>1 keV) have 
a high penetration factor. Th ey can pierce through colu mns of 
gas of up to N ^ 10^** cm - (iMeiierink & Spaansll2005h . Since 
the X-ray source is in the center of the simulation box, material 
that lies behind high density gas, with column densities beyond 
10"^ cm"^, can be strongly shielded from the radiation. This ac- 
tually occurs in our simulations. Our initiahzed model cloud has 
a maximum column of N = 2.0 ■ 10^^ cm~^ (p = 10^ cm""* and 
2rcioud = 2.0 • 10"^ cm). Some of the phase diagrams show a 
secondary line at low densities that exhibits a sharp drop to low 
temperatures because of the cold shielded gas. In each of the 
Figs. |6] to [17] the X-rays increase from top to bottom. 



From the phase diagrams, it is immediately clear that the 
gas temperatures decrease with increasing density and that the 
temperatures are higher for the higher X-ray fluxes at any point 
in density. The latter is a direct consequence of the efficient 
Coulomb heating by X-rays. X-rays do not only heat the system, 
but also allow new paths for cooling to proceed along. X-rays in- 
crease the ionization fraction of the species and new molecules 
are formed. This results in a higher cooling efficiency with in- 
creasing density, hence, the decreasing temperature trend. We 
can also see a spread in the temperature as a direct effect from 
the different column densities throughout the cloud. The spread 
seems to be larger for the higher X-ray fluxes, however, this is 
not resulting from a larger column density range, but is merely 
due to the larger range in temperatures from the increased X- 
ray flux. An interesting feature is that there is a range in densi- 
ties around n ~ 10^ cm"^ where the temperature lingers around 
a few 100 K. The cooling process is slowed down here due to 
LTE effects (thermalization and line trapping) until the densities 
are high enough so that cooling can start to be effective again 
through gas-dust coupling. 

From the color in the phase diagrams we can infer that most 
of the gas mass lies at high densities, but for higher X-ray fluxes 
a lot of mass also lies in the low density, high temperature 
regime. This is the area where the cloud is directly irradiated by 
X-rays and the cloud evaporates. As a consequence, the cloud 
size and its mass are reduced by this. This has repercussions for 
the final star-formation efficiency. When comparing an isother- 
mal model against the model with the highest X-ray flux, for 
^CR = 1, we see that the SFE is reduced from 9.4% to 8.2% 
for a 10^ Mo black hole at t = 3tff, and from 12.3% to 7.0% 
for a 10^ Mq black hole (Figs. [TOl and [TTTi. The same behavior 
is found for higher cosmic ray rates (discussed in section |42]i. 
Despite some fluctuations due to the randomness of the system, 
it is interesting to see that the efficiency generally increases a lit- 
tle with our lowest X-ray flux compared to the isothermal mod- 
els. It seems that a mild X-ray flux of 5. 1 erg s"' cm"- or less ac- 
tually enhances the star-formation efficiency. See Table|2]for the 
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list of results. Although gas temperatures are higher and the SFE 
decreases when going to higher X-ray fluxes, the star-formation 
rate, SFR = dSFE/dt, is still high early on. We can see from the 
SFE plots that star formation is initiated at about the same time, 
tonset - 0.6tff, irrespective of X-rays. We find that the quench- 
ing of star formation due to X-ray heating is balanced by the 
increased densities resulting from the thermal pressure that the 
radiation field creates. An ionizing pressure front compresses the 
cloud which leads to very efficient star formation, similar to, if 
not more than, the colder, isothermal, cloud models. Fig.[T]shows 
a 2D slice to this effect. 




22e+04 



4.7e+02 



Fig. 1. Density morphology of model M19, i.e., with Fx - 
160 erg s"' cm"-, at t=2/3tff. The image shows a slice through 
the XY-plane of the center of the cloud. The color represents the 
number density (cm^-^) and the axes are in parsec. The arrows 
represent the direction of radiation emanating from the black 
hole, which is located at the upper left side. 



We see that the IMFs tend to get flatter for the higher X- 
ray fluxes. The best power-law fit gives us a nice Salpeter slope 
with Fpowfit = -1.34 for model MOl, i.e., without X-rays, and 
a slope of Fpowfit = -1.31,-1.24,-1.01 for models M07, M13, 
and M19, i.e., for fixed parameters but with increasing X-ray 
flux. A lognormal distribution also seems to fit most of the 
IMFs quite well. The same trend is visible for all other cloud 
conditions. This implies that the IMF shape depends on how 
much the cloud is irradiated. Another effect is that the turn- 
over mass shifts toward higher masses with increasing flux, from 
Mturn = 0.54 Mq, at the lowest X-ray flux (M07), to about 0.70 
and 1.51 Mq for the higher fluxes (M13 and M19). However, the 
corresponding isothermal fiducial model (MOl) seems to have a 
bit higher turn-over mass, 0.66 Mq, than our lowest X-ray flux 
(5.1 erg s"' cm"^) model. Again, a similar trend is seen for all 
cloud models. 



4.2. Effects of cosmic rays 

Cosmic rays can easily penetrate a molecular cloud and heat 
the gas uniformally. As such, the y set the minimum attain- 
able temperature in these systems dGoldsmith & Langeijll978l: 
|i3_ergin & Tafalla 2007). For our model cloud conditions, which 
have cosmic ray rates of 1, 100, and 3000xGalactic (^cR.Cai = 
5 X 10"'^ s"'), the minimum temperatures are 10 K, 50 K, and 



200 K (iPapadopoulos et"al ]|201 lUMeiierink et al.ll201 ih . In each 
of the Figs. |6] to [IT] the cosmic ray rate increases from left to 
right. 

In contrast to the X-ray heating, cosmic rays clearly delay 
the onset of star formation. The first sink particles are formed 
much later with increasing ^cr, see Fig. |2] The number of sink 
particles formed is also drastically reduced by up to one order 
of magnitude, see Table |2] Despite this, once the cloud starts 
forming stars, it can do so rapidly and massively. The much 
higher mass each star obtains compensates for the loss in the 
number of stars. In the end, the SFEs between the different ^cr 
runs are not much different. This can be explained by the fact 
that the Jeans mass scales strongly with the temperature and that 
the cloud fragments into star-forming cores that are more mas- 
sive when the ambient temperature is higher The stars that form 
out of them are therefore also heavier. These massive stars ac- 
crete more matter, since there is less competition around, and 
they tend to grow more. This result is also in agreement with the 
idea about fragmen tation-induced starvation (iPeters et al.ll2010l: 
iKlessen et al l20Toh . But, we have a lack of fragmentation in this 
case and hence very little starvation. 



1.2 



1.0 - 



-0.8h 

C/3 

H — 

O 

o 0-6 

C/3 

O 

0.4 



• CcR=1xGal. 
oCcR=100xGal. 
A ?cR=3000xGal. 








1 10 
X-ray flux (erg s'^ cm'^] 



100 



Fig. 2. Onset of star formation for the runs with a 10^ M© black 
hole. In this figure, the star-formation initiation time is plotted 
against X-ray flux for three different cosmic ray rates. Higher 
cosmic ray rates delay the onset of star formation, while increas- 
ing the X-ray flux counteracts this. Note that the points on the 
left side are for the isothermal models which have an X-ray flux 
ofO. 



In the absence of X-rays, one obvious result is that high cos- 
mic ray rates strongly inhibit the formation of sub-solar mass 
sink particles, which follows from the aforementioned Jeans 
mass argument. The smallest mass we obtain in these runs is 
ab out 1 solar mass. This is i n agreement with, and predicted 
bv. IPapadopoulos et al.l (1201 ih . From the IMF plots, we see that 
the effect of cosmic rays on the slope of the mass function is 
minimal. The lognormal shape is also not much affected, al- 
though, given the relative low number of points at higher ^cr 
values, one cannot firmly conclude this. It is thus so that there is 
a general shift in the mass distribution toward higher masses. 
The turn-over mass scales toward higher masses as well and 
has shifted from 0.66 Mq (MOl) to 2.37 Mo (M02) for the 
Mbh = 10^ Ma runs, and from 0.47 Mq (M04) to 2.60 Mq (M05) 
for the Mbh - lO*" Mq runs. This is an increase of a factor of 
3.6 and 5.5, respectively. Between these runs, the Jeans mass in- 
creases by a factor of 5-^^^ = 11.2 (10 K to 50 K) at the same 
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Table 2. List of results for each model 



Model 


IMF slope 


Turn-over mass 
[Mo] 


Onset of SF 
[tff = 10^ yr] 


Final SFE 

[Mp,n7M,„,.] 


Accreted fraction 


Number of particles 


MOl 


-1.34 


0.66 


0.47 


0.094 


0.606 


247 


M02 


-1.44 


2.37 


0.87 


0.092 


0.691 


110 


M03 


— " 


4.25" 


1.16 


0.089 


0.747 


27 


M04 


-1.30 


0.47 


0.35 


0.123 


0.801 


262 


M05 


-1.24 


2.60 


0.67 


0.125 


0.795 


63 


M06 


— " 


5.95" 


0.99 


0.117 


0.725 


21 


MOV 


-1.31 


0.54 


0.43 


0.099 


0.609 


350 


MOB 


-1.08 


0.65 


0.50 


0.095 


0.693 


227 


M09 


-1.07 


0.65 


0.60 


0.077 


0.653 


211 


MIO 


-1.04 


0.43 


0.36 


0.125 


0.768 


276 


Mil 


-1.05 


0.52 


0.46 


0.117 


0.708 


137 


M12 


-0.78 


0.26 


0.52 


0.116 


0.762 


104 


M13 


-1.24 


0.70 


0.52 


0.102 


0.607 


286 


M14 


-1.02 


0.73 


0.59 


0.092 


0.597 


211 


M15 


-0.92 


0.77 


0.67 


0.093 


0.619 


171 


M16 


-1.16 


0.68 


0.42 


0.104 


0.598 


202 


M17 


-1.01 


0.65 


0.52 


0.100 


0.545 


153 


M18 


-0.58 


0.52 


0.59 


0.094 


0.592 


104 


M19 


-1.01 


1.51 


0.65 


0.082 


0.637 


138 


M20 


-0.91 


1.62 


0.65 


0.076 


0.620 


134 


M21 


-0.85 


1.43 


0.67 


0.073 


0.621 


111 


M22 


-0.82 


1.12 


0.55 


0.070 


0.647 


80 


M23 


-0.70 


1.91 


0.57 


0.071 


0.623 


79 


M24 


-0.62 


1.36 


0.62 


0.073 


0.627 


58 


M25 


-1.32 


1.20 


0.74 


0.014 


0.476 


84 


M26 


— 


— 


0.80 


0.001 


0.681 


3 


M27 


— * 


— * 


— ' 


— * 


— * 





M28 


-1.30 


1.62 


0.66 


0.037 


0.637 


108 


M29 


-1.28 


1.86 


0.66 


0.024 


0.697 


59 


M30 


-1.07 


1.53 


0.65 


0.013 


0.618 


49 


M31 


-1.48 


0.55 


0.49 


0.103 


0.638 


262 


M32 


-1.31 


2.30 


0.90 


0.099 


0.728 


106 


M33 


— " 


4.30° 


1.16 


0.093 


0.850 


19 


M34 


-1.23 


0.43 


0.40 


0.101 


0.568 


332 


M3j 




U.D / 


U.j4 


U.UOJ 


U.DDS 




M36 


-1.06 


0.65 


0.60 


0.096 


0.638 


198 


M37 


-1.05 


0.59 


0.51 


0.107 


0.650 


281 


M38 


-1.04 


0.66 


0.56 


0.101 


0.647 


207 


M39 


-0.99 


0.88 


0.67 


0.101 


0.748 


163 


M40 


-1.40 


1.58 


0.66 


0.084 


0.641 


142 


M41 


-1.06 


1.47 


0.66 


0.082 


0.648 


132 


M42 


-0.79 


1.52 


0.67 


0.069 


0.614 


103 



" Too few particles have formed here to properly construct an IMF. 
* No particles have formed in this simulation. 



densities. So, for ^cr = 1 to 100 times Galactic, it seems that 
the turn-over mass scales with roughly the temperature rather 
than the Jeans mass. We find that very massive sink particles are 
formed in the highest cosmic ray simulations. The most massive 
sink particle that formed in model M03 has a mass of 350 solar 
masses. 

4.3. Effects of UV 

In PDRs, UV radiation heats the gas up to a few thousand K 
and will dominate the chemistry of cloud surfaces. UV radia- 
tion can photodissociate molecules, which will lead to the de- 
struction of efficient coolants, and heats through photo-electric 
emission from (small) dust grains. UV radiation is strongly at- 
tenuated above columns of N = 10~^ cm"^, at solar metallic- 
ity. In our simulations, the model cloud has a starting column of 
N - lO^-'cm"- (center to edge). This grows to N - 1x10^"^ cm"^, 
and beyond, when collapse is initiated and the first stars are 



about to form. Hence, it is not expected that UV radiation will 
have a direct effect on the results. Only by heating the low den- 
sity regions, n < lO'* cm"-', and by pressurizing the medium 
where the molecular cloud is embedded in, the results might be 
influenced. We assume that there is a uniform background UV 
field. The results for the model clouds in UV dominated envi- 
ronments are shown in Figs. |9] US andfTTl 

We have chosen to simulate the models with increased UV 
radiation for a single black hole mass, Mbh = 10^ Mq, as can 
be seen from the table (Table [T]i. We have done so because our 
initial studies immediately showed that the direct effect of an ex- 
ternal U V radiation field on our model clouds is modest. Follow- 
up simulations confirmed that the column density (>10^~ cm"^) 
of the molecular cloud is, as predicted, already too high for UV 
radiation to penetrate even the edge of the cloud (<0.01 pc) and 
that turbulent eff'ects and other instabilities do not provide the 
opportunity for UV radiation to do much damage. When X-rays 
are present, the gas heating is dominated by them. This is a con- 
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Fig. 3. Comparison of IMFs of the UV runs against their non-UV counterparts. The left panel shows the IMFs of the models MOl 
and M31 (no X-rays). The right panel shows the IMFs of the models M19 and M40 (highest X-ray flux). For comparison purposes, 
the Salpeter IMF (green dashed line) and the Chabrier IMF (blue dot-dashed line) are shown, as fitted to our fiducial case. 



sequence of the much higher (10-50%) ionization heating effi- 
ciency compared to photo-electric heating (0.1-1%) by UV ir- 
radiated dust grains. However, UV radiation does heat the low 
density gas, therefore, it influences the late time accretion. 

The phase diagrams for the UV runs confirm that the radia- 
tion does not penetrate the natal cloud beyond a few magnitudes 
of visual extinction at one free-fall time. At low densities, we see 
that UV heats the gas to several hundred and upto a few thousand 
K. Above 10^ cm"^ it has no impact. In the models together with 
X-rays, even with the lowest X-ray flux Fx = 5.1 erg s"' cm"^, 
there is little difference in the phase diagrams, as compared to the 
non-UV runs. Only the low density < 10"^ cm"^, unshielded (in 
the UV) regions are heated up by UV. We also see that the SFE 
plots are not much different than their UV deficient counterparts. 
The slope of the IMF, on the other hand, is more sensitive to vari- 
ations. We see that massive stars accrete less matter during the 
late stages of the run, i.e., t > 2tff, as the temperature is higher 
at low densities. This while most of the late time accretion nor- 
mally comes from the low density, cold regions. Especially for 
tiie highest Fx runs, the slope of the IMF becomes slightly more 
steep due to a decrease in the mass of the most massive stars, in 
the presence of a UV background. In Fig. [3] we plot the IMFs 
for four model runs, MOl versus M3 1 and M19 versus M40, that 
is, with and without UV, together in the same figure to highlight 
this eff'ect. 

4.4. Effects of shear 

In the presence of strong gravity, the gravitational pull on each 
side of a bound extended object will cause a differential acceler- 
ation. If the kinetic energy resulting from the velocity difference 
is close to or on the order of the gravitational binding energy of 
the object, it will play an important role in the dynamics of the 
cloud. The object will be torn apart if the corresponding kinetic 
energy is much higher than its binding energy. 
Our model cloud has a binding energy of 

3GM-, , 

Eb,nd = = 1.02 X lO'*^ erg, (10) 

-^I'cloud 

where Mdoud and rdoud are the cloud mass and radius. The max- 
imum shear for the models with the three different black hole 



masses, for a cloud that is at 10 pc distance, results in kinetic 
energies of Eshear = [4 x 10"*^ 4 x 10'**', 4 x 10^^] erg, with 
AVshear = [0.7,2.2,7.0] km/s. The velocity difference result- 
ing from the gravitational stresses of the two higher black hole 
masses, Mbh - [10^, 10^] M©, is on the order of the applied ini- 
tial turbulence (5 km/s), and will therefore play an important role 
in the evolution of the cloud. Whereas the velocity shear due to 
a black hole with Mbh = 10^ M© will have very little impact 
on the dynamics of the simulation. Figs.|6][T0l and [14] display 
the models with Mbh = 10^ Mq. Figs. 01111 and [15] display the 
models with Mbh = 10'' Mq. While, Figs.[8][l2] and[T6]show the 
models with Mbh - 10'^ Mq. 

From the phase diagrams we infer that the column density 
range is more extended for the lower shear runs. A salient fea- 
ture is the dual phase structure that can be seen in the phase 
diagrams at low densities. This is because the cloud retains its 
size and shape for a longer period of time when there is less 
shear, developing a large column density range for the same den- 
sities. Moreover, we see from the density images that the cloud 
becomes more asymmetric with increasing Mbh. The cloud is 
stretched in the direction of the orbit but compressed perpen- 
dicular to it, see Fig. ]4] Due to their morphological effect on the 
cloud, the gravitational stresses are the cause for the reduced col- 
umn density range. This causes that in the presence of X-rays, a 
larger surface is directly irradiated by radiation with less attenu- 
ation throughout the cloud. 

We find that the final star-formation efficiencies are reduced 
for higher black hole masses. Furthermore, we see that the star- 
formation rate has dropped as well, while the onset of star for- 
mation is delayed with increasing shear. However, when we look 
at the numbers of sink particles formed, we find that they are 
remarkably higher for large shear, that is, comparing Mbh = 
10^ Mq against 10^ Mq. The simulations with the highest shear, 
Mbh = 10*^ Mo, strongly inhibit the formation of all stars. See 
Table [2] for details. This increase in the number of sink parti- 
cles, while the SFE does not, shows that the formation of low- 
mass stars is favored for the case with Mbh = 10^ Mq. When we 
look at the IMF plots, we do see that for the higher shear runs, 
the mass function generally comprises more low-mass stars and 
fewer high-mass stars . This result is in agreement with grav otur- 
bulent fragmentation dKlessen & Ballesteros-ParedeJ2004l) . 
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Fig. 4. Density slices through the center of the cloud at t = 2/3tff for three different black hole mass models. The clouds are stretched 
in the direction of motion with increasing black hole mass, but compressed in the direction perpendicular to this. The black hole 
masses are: 10^ Mq for model M04 (left panel), 10^ Mq for model MOl (middle panel), and 10^ Mq for model M25 (right panel). 
The white arrow shows the direction of the orbital motion. 



4.5. Accretion onto sink particles 

Generally most of the mass growth of the sink particles comes 
from accretion. Proto-stars gain about 66% of their mass in this 
way. This gives us an average accretion rate of ~10"^ Mo/yr. 
Due to the restricted time resolution of our checkpoint files, im- 
mediate merging of newly created sink particles within 1/30* 
of a free-fall time (~3xlO^ yr), is considered as accretion. 
However, this should not occur very often, since only about 15% 
of the particles merge during a run and this is spread throughout 
the whole simulation. The fraction of accreted mass at t = 3tff 
is, on average, facer - 0.683, when there is little black hole shear 
(Mbh - 10*' Mq). This reduces to facer - 0.642 when we increase 
Mbh to 10^ Mq. Similarly, when we increase the X-ray flux from 
0, 5.1, 28 to 160, the accreted mass fraction drops from 0.774, 
0.746, 0.579 to 0.632 for Mbh = 10'' Mq and from 0.681, 0.652, 
0.608 to 0.626 for Mbh = 10^ Mq. See Table |2]for a complete 
list of accreted mass fractions. Note that the highest X-ray flux 
has a bit higher mass accretion fraction than the second high- 
est flux. This may be a consequence of the interplay between 
reduced accretion due to the increase in temperature, as sound 
speed scales with T'^-, and the increase in accretion due to the 
larger gas reservoir available because of the lower star-formation 
efliciency. Apparently, the balance tips towards more accretion at 
the highest X-ray flux. The fractions for the runs with the highest 
shear are too volatile to consider, since very little star formation 
has occurred in these runs. 

From Eq. |5] we see that the accretion rate scales inversely 
with the sound speed and gas velocity to the third power This 
is a strong scaling and is the dominant factor for the decrease in 
accretion as gas turbulence increases. The strong differential ve- 
locities due to shear from the black hole cascade down to smaller 
scales and yield turbulence. The turbulence still decays in time, 
but the decay is less rapid compared to when there is less shear 
Thermal pressure as a consequence of the increased tempera- 
tures for the X-ray models also contributes to higher velocities 
and decreases the accretion rate. We do not see a clear correla- 
tion between accretion and cosmic rays, though. There can be 
two reasons for this. First of all, there is no gradient in pressure 
throughout the cloud due to the increase in cosmic ray ionization 



rate (they heat all gas uniformly). Second, the heating by cosmic 
rays is relatively modest, yielding much lower temperatures than 
that of UV or X-ray heated gas. 

For the UV runs, we find notable diff'erences after two free- 
fall times. Due to UV heating, accretion onto sink particles from 
the low density environment, n < 10'* cm"^, is drastically re- 
duced. This is normally the main source for particles to gain 
mass at the later stages of the simulation. The more massive 
sink particles are also more strongly affected. Consequently, the 
models that produce a top-heavy IMF, these are the runs with 
increased X-ray flux, are restrained. A slightly more steep IMF 
slope is the result. The turn-over mass and the SFEs, on the other 
hand, are not significantly affected. 



5. Conclusions and discussion 

By means of numerical simulations, we have performed an ex- 
tensive study on the effects of various kinds of irradiation (X- 
rays, cosmic rays, and UV) and shear on collapsing molecular 
clouds in active galactic environments. We have analyzed how 
the star-formation efficiencies and the initial mass functions are 
affected for 42 different cloud models. Our general result is that 
the IMF deviates from a Salpeter shape in extreme environments. 
Still, for a mixture of ambient conditions, their effects on the 
mass function somewhat balances out and a deviation from a 
Salpeter shape is less pronounced. We evaluate the differences 
and similarities in detail. All the quantifiable results are listed in 
Tabled 

The IMF nicely follows a lognormal distribution with a 
Salpeter slope for the isothermal runs. Our fiducial model MOl 
and model M04, which have similar conditions to those of the 
Milky Way, are in good agreement with a Chabrier IMF and 
match a Salpeter slope excellently, as can be seen in Fig. [14] 
Increasing the cosmic ray rate or the black hole shear does not 
change the slope of the mass function by much, except perhaps 
for the highest ^cr and the highest Mbh. Due to the insufficient 
number of sink particles, it is not possible to make a good fit for 
these cases. 
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The interplay between X-ray, shear and cosmic rays is as 
follows. The power-law slope of the IMF flattens and becomes 
non-Salpeter when we increase the X-ray flux from to 160 
erg s"' cm"^. We find that the slope drops from around Fpowfit - 
-1.35 (Fx = 0) to about Fpowfit = -1 (Fx = 160) for the runs 
with Mbh = 10^ Mo, and to about Fpowfit = -0.8 (Fx = 160) 
when Mbh = 10*" M©. We see that the flattening is less promi- 
nent when the shear is stronger We attribute this effect to the 
change in shape of the cloud, see Fig. |4] When gravitational 
shear is stronger, the irradiated face of the cloud is larger due 
to stretching of the cloud in the direction of the orbit. A bigger 
region is thus compressed by thermal pressure. We also see that 
the cloud obtains a rotational motion. Its rotational time is com- 
parable to, though slightly faster than, the orbital time. This rota- 
tion is slow, however, it does cause the regions that were initially 
hot, to become shielded from X-rays. These regions, which had 
their densities raised by thermal compression now cool down 
quickly in the shade and form low-mass stars easily. This while 
the model clouds in the lowest Mbh - 10^ M© runs face the 
radiation source on the same side for a much longer period, as 
their orbital time is about three times larger. This argument is 
supported by the behavior of the models with the highest shear. 
For these, the isothermal cosmic ray runs do not form sink par- 
ticles at all, however, in the presence of X-rays, star formation 
does still occur thanks to compression. Cosmic rays, on the other 
hand, enhance the flattening of the power-law slope. At the high- 
est ^CR. the IMF slope drops further from about Fpowfit =-1 to -0.8 
for the Mbh = 10^ M© runs, and from Fpowfit =-0.8 to -0.6 for the 
Mbh - 10^ Mq runs. We can understand this behavior since there 
are regions in the cloud that are shielded from X-rays and where 
the gas behaves in an isothermal manner. These places will form 
more massive stars because of the higher temperatures. Since 
we plot all the stars together in the IMF plot, a steep overall IMF 
slope is created. 

We find that the turn-over mass shifts towards higher masses 
for increasing X-ray fluxes and cosmic ray rates. A simple Jeans 
mass argument is sufficient to explain this behavior. Since the 
Jeans mass strongly depends on the temperature (Eq. O, and 
this sets the minimum mass for the fragmentation scale, stars 
that form in massive cores tend to be more massive as well. The 
effect of cosmic rays is less prominent in the presence of X-rays 
due to the softer equation of state, as mentioned earlier. This, in 
the end, makes the turn-over mass less strongly dependent on the 
Jeans mass, and in part explains the constancy of the character- 
istic mass in many star-forming regions (Elmegreenet al. 2008^ 
lBatgi2009l) . 

X-rays reduce the final star-formation efficiency by up 
to 40% in the presence of a strong X-ray flux. Fx - 
160 erg s"' cm"^. This effect is very modest for the lowest X- 
ray flux. Fx =5.1 erg s"' cm"^. In fact, there is a slight increase 
in the SEE for the lower cosmic ray rates. For X-rays to have any 
significant impact on to the efficiency of star formation, a flux of 
at least Fx > 5.1 erg s"' cm"^ is needed. The main reason for 
the reduced efficiency is that the cloud evaporates from the irra- 
diated side, thereby reducing the mass of the molecular cloud. 
This, however, also increases the densities at the same side due 
to an ionizing compression front such that the star-formation rate 
remains high. 

We also see a reduction in formation efficiency when we in- 
crease the black hole shear We do this by increasing the black 
hole mass while keeping the distance of the cloud fixed at 10 pc, 
thereby changing the Mbh/dj^j^ ratio. For the runs with the high- 
est shear, Mbh = 10^ Mq, star formation is almost completely 



quenched. In this case, the kinetic energy from the velocity di- 
vergence that the black hole injects is greater than the binding 
energy of the cloud. However, we find that stars can still form 
when thermal compression by X-rays enhances cloud collapse, 
albeit, with very low efficiency. The effect of shear on star for- 
mation is also more dramatic than that of the X-rays. Besides 
reducing the efficiency, the rate at which stars form is also af- 
fected. From the shape of the curves of the SEE plots in Figs. 
[T0]to[T3] we see that the increase in SEE in time is less steep 
and somewhat more irregular for higher Mbh. The cause for this 
effect is the enhancement of the turbulence from the shearing 
motion cascading down to smaller scales. We see that the turbu- 
lence remains strong, FWHM = 4-5 km/s on the scale of the 
cloud at one free-fall time, for the Mbh = 10^ and 10^ Mq runs, 
while the turbulence decays to about 1 km/s for Mbh = 10^ Mq. 
In addition to this, stars start to form later with increasing shear 

Cosmic rays also delay the onset of star formation. 
Interestingly, the star-formation efficiency and the star-formation 
rate are not affected. However, the number of sink particles is 
drastically reduced with increasing ^cr- The average mass per 
sink particle, on the other hand, rises comparably, while the for- 
mation of low-mass stars is strongly inhibited. There is a signifi- 
cant side-effect from X-rays that counteracts this. When we look 
at the difference between the isothermal runs and the X-ray runs, 
for any ^cr, we see that many more sink particles are formed 
when there is an X-ray source. Especially the number of low- 
mass stars is higher for the X-ray runs. This can be inferred form 
the shape of the mass function in Figs.[T4]to[T7] One would nor- 
mally expect that due to the increased gas temperatures from X- 
ray heating, for the same cosmic ray rate, star formation would 
be strongly suppressed. However, the equation of state plays a 
crucial role here. Despite the higher temperatures, the compress- 
ibility of the gas, y being less than unity, causes the molecular 
cloud to fragment more easily. This while the cosmic ray runs 
enjoy an isothermal, y = 1, equation of state. Therefore, low- 
mass stars are able to form with much less effort. The delay in 
star formation due to cosmic ray heating is also counterbalanced 
in the presence of X-rays. This effect is more prominent with 
increasing ^cr, see Fig.|2l A time difference in the onset of star 
formation of more than one-third free-fall time (> 3.3 x 10"* s) 
can be seen for models with ^cr = 3000 x Galactic. 

For our models with UV, we find that it has modest impact 
on the results. This was expected, since our model cloud has a 
column of N ^ lO^"' cm"^ at the start of the simulation, and in- 
creases its (column) density with time. UV radiation is strongly 
attenuated above columns of N = 10'^ cm"^, for solar metallic- 
ity. We presented the models with UV for Mbh =10^ Mq and 
Fuv = 10^^ Go only, to asses whether increased external pres- 
sure and the higher cloud edge gas temperatures can play an im- 
portant role in the evolution of the cloud. We find that by increas- 
ing the temperature for the lower densities, n < 10"* cm"^, ac- 
cretion onto particles is reduced at later times, i.e., t > 2tfl^. This 
impacts the growth of the massive particles later on. As such, the 
slope of the IMF, especially in the presence of X-rays, is some- 
what more steep, in general. However, we find that the change 
in the SEEs are almost negligible when we have an isotropic UV 
field with a strength of 10~^ Go, as compared to the non-UV 
runs. For cloud conditions such as in this work, we find that UV 
radiation can still be of importance, through secondary channels, 
to the formation of stars. 

Accretion rates onto point particles in a turbulent medium 
may be different than the rates (Eq. |2ll used in this work as 
iKrumholz et al.l (l2006h point out. Our medium is close to homo- 
geneous on scales of a few parsec. On the other hand, outflows 
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and jets from young stars could reduce the net accretion rates on 
smaller scales. Therefore, using a homogeneous flow approxi- 
mation might overestimate the accreti on rates and we c o nsider 
them upper limits. In addition to this, IStamatellos et al.l (1201 ih 
show that episodic accretion for low-mass stars with bursts of 
radiative feedback can also aff'ect the fragmentation properties. 
We have not taken such effects into consideration. We compare, 
however, our model results against each other In this way, our 
findings are less sensitive to the above mentioned issues, and 
relative changes in the IMF are well defined. 

The velocity shear through the cloud due to black hole grav- 
ity decreases as the cloud contracts. With time, the model cloud 
in our simulations shrinks in the direction perpendicular to the 
orbit due to self-gravity and turbulent motions. On average, the 
model cloud, shrinks by a factor of two, at most, in one free-fall 
time. This reduces the shear by the same amount as Av scales 
with rcioud- Thus, the injected energy decreases with time, but 
not strongly. 

In this study, the maximum resolution allowed was 4096^ 
cells. This in itself is quite a good resolution, and to follow it for 
long dynamical timescales, like we do, makes it difficult to go 
much higher Since we use an adaptive grid, it is more impor- 
tant to properly resolve the grid dynamically. A proper means 
to do so is to re solve the Jeans l ength adequately everywhere in 
time and space. FlYuelove et"an (Il997h say that the Jeans length 
should be resolved by at least 4 cells, but in the case of (M)HD 
simulations in magnetic fields, a minimum reso lution of 30 cells 
is required according to lFederrath et al.l (1201 ih . We resolved the 
Jeans length in our simulations by at least 10 cells. This actually 
means that the minimum resolution always lies somewhere be- 
tween 10-20 cells. However, our resolution criterion is so strict 
that we resolve a whole block of 512 cells once even a single 
cell comes close to the minimum of 10 cells, while we only de- 
refine if all the cells within the block have stretched beyond 25 
times the Jeans length. This in effect makes the general resolu- 
tion in our runs much higher than the 10 cells that we mention, 
and approximately about 50-60 cells on average. 

Stellar multiplicity is not covered in this work and binaries 
do not form naturally from the sink particle routine. In order to 
achieve that, cloud collapse should be followed to the highest 
grid resolutions in order to resolve the gas dynamics properly 
and allow proto-stars to form. This is very difficult to perform 
numerically. If we were to consider the sink particles in our work 
as binaries, then the found shapes of the IMFs would not be dif- 
ferent, and would only be rescaled to lower masses. As a test, we 
have re-run 2 of our simulations without sink particle merging. 
We found that the merging routine does not influence the results 
much. Fig.|5]shows a plot of model MOl without merging of sink 
particles. 

In the future, we intend to perform studies for conditions 
much closer to the black hole, i.e., within 1 pc of galactic nu- 
clei, to model massive stars as observed in Sgr A* and in M3 1 . 
We also plan to incorporate our results into simulations of AGN 
evolution and intend to look into the stability of the entire disk. 
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Fig. 5. The IMF of model MOl without sink particle merging. In 
this image, the Salpeter IMF (green dashed) and the Chabrier 
IMF (blue dot-dashed) are displayed as fitted to our fiducial 
model. A linear fit and a lognormal fit are shown as purple and 
red lines. 
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Fig. 6. Phase diagrams for the models with Mbh = 10^ and UV=0. Temperature is plotted against number density at t = tff = 
10^ yr. The diagrams are gridded into 75~ cells with the weighted masses of the points depicted in color. Red represents a mass of 
10 Mq or above and blue is for 10"^ M©. Isothermal conditions yield a flat profile. From left to right, the cosmic ray rate increases 
from 1, 100 to 3000 x Galactic. From top to bottom, the X-ray flux increases from 0, 5.1, 28 to 160 erg s"' cm"^. 
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Fig. 7. Phase diagrams for the models with Mbh = lO*" and UV=0. Temperature is plotted against number density at t = tff = 
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Fig. 10. Star-formation efficiencies for the models with Mbh - 10^ and UV=0. The ratio of the total sink particle mass over the 
total initial gas mass is plotted against time (in free-fall units). From left to right, the cosmic ray rate increases from 1, 100 to 3000 
X Galactic. From top to bottom, the X-ray flux increases from 0, 5.1, 28 to 160 erg s"' cm"^. The total number of sink particles 
formed during the run is given in the upper left corner of each panel. 
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Fig. 11. Star-formation efficiencies for the models with Mbh - 10'' Mq and UV=0. The ratio of the total sink particle mass over the 
total initial gas mass is plotted against time (in free-fall units). From left to right, the cosmic ray rate increases from 1, 100 to 3000 
X Galactic. From top to bottom, the X-ray flux increases from 0, 5.1, 28 to 160 erg s"' cm"^. The total number of sink particles 
formed during the run is given in the upper left corner of each panel. 
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Fig. 12. Star-formation efficiencies for the models with Mbh = 10^ Mq and UV=0. The ratio of the total sink particle mass over 
the total initial gas mass is plotted against time (in free-fall units). Note that the y-axis range in this figure differs from the other 
(SFE) figures. From left to right, the cosmic ray rate increases from 1, 100 to 3000 x Galactic. From top to bottom, the X-ray flux 
increases from to 160 erg s"' cm"^. The total number of sink particles formed during the run is given in the upper left corner of 
each panel. 
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Fig. 13. Star-formation efficiencies for the models with UV. The UV flux used in these models is 10^-^ Go- Similar to the figures 
before, Figs.[T0]to[T2l the images show the ratio of the total sink particle mass over the total initial gas mass and is plotted against 
time (in free-fall units). From left to right, the cosmic ray rate increases from 1, 100 to 3000 x Galactic. From top to bottom, the 
X-ray flux increases from 0, 5.1, 28 to 160 erg s"' cm"^. These UV models are simulated for a a black hole mass of Mbh - 10^ Mq. 
The total number of sink particles formed during the run is given in the upper left corner of each panel. 
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Fig. 14. Initial mass functions for the models with Mbh = 10^ and UV=0. The images display the time averaged IMFs between 
1 and 3 free-fall times, where tff - 10^ yr. From left to right, the cosmic ray rate increases from 1, 100 to 3000 x Galactic. From top 
to bottom, the X-ray flux increases from 0, 5.1, 28 to 160 erg s"' cm"-. In each image, for comparison purposes, the Saipeter IMF 
(green dashed) and the Chabrier IMF (blue dot-dashed) are displayed as fitted to our fiducial model. Two best fits are applied to the 
data, a linear fit and a lognormal fit, and are shown as purple and red lines. With the exception of the lognormal fit, the slopes above 
the turn-over mass are given in the upper left corner. 
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Fig. 15. Initial mass functions for the models with Mbh = 10^ and UV=0. The images display the time averaged IMFs between 
1 and 3 free-fall times, where tff - 10^ yr. From left to right, the cosmic ray rate increases from 1, 100 to 3000 x Galactic. From top 
to bottom, the X-ray flux increases from 0, 5.1, 28 to 160 erg s"' cm"-. In each image, for comparison purposes, the Saipeter IMF 
(green dashed) and the Chabrier IMF (blue dot-dashed) are displayed as fitted to our fiducial model. Two best fits are applied to the 
data, a linear fit and a lognormal fit, and are shown as purple and red lines. With the exception of the lognormal fit, the slopes above 
the turn-over mass are given in the upper left corner. 



24 



100 



M25 



S. Hocuk & M. Spaans: Star formation in AGN 
M26 



■° 10 



0.1 



• \ • ■ ■ ■ 


■ ' ' 1 ' 


— 1 — 1 1 ; 


:rsa, =-i.3&\ 

. Tcha =-1.30 




Power-law fit ', 

Lognormal fit ■ 

Salpeter IMF ■ 

CliabrierlI\^F' 



















100 



10 



1.0 

M(Me) 



10.0 



r„ - .1.30 



0.1 



Power-law fit 
Lognormal fit 
Salpeter IIVIF 
Ctiabrler IIVIF 



100 



10 



1.0 

M(Me) 



10.0 



M27 



1 1_ 
0.1 



= -1 .35 
= -1 .30 



■ Power-law fit 

■ Lognormal fit 
Salpeter IIVIF 
Ctiabrler IIVIF 



1.0 

M(Me) 



10.0 



100 



10 



0.1 



M28 



M29 



M30 



= -1.30 
= -1.35 
= -1.30 



. Power-law fit 

■ Lognormal fit 
Salpeter IMF 

■ Chabrier IMF 




100 



10 



= -1 .28 
= -1 .35 
= -1 .30 



. Power-law fit 

■ Lognormal fit 
SaTpeter IMF 

■ Ctiabrier IMF 




100 



10 



= -1 .07 
= -1 .35 
= -1 .30 



. Power-law fit 
■ Lognormal fit 
Salpeter IMF 
. CfiabrierlMF 




1.0 

M(Me) 



10.0 



0.1 



1.0 

M(Me) 



10.0 



0.1 



1.0 

M(Me) 



10.0 



Fig. 16. Initial mass functions for the models with Mbh = 10^ M© and UV=0. The images display the time averaged IMFs between 
1 and 3 free-fall times, where t^ = 10^ yr. From left to right, the cosmic ray rate increases from 1, 100 to 3000 x Galactic. From top 



to bottom, the X-ray flux increases from to 160 erg s" cm" . In each image, for comparison purposes, the Salpeter IMF (green 
dashed) and the Chabrier IMF (blue dot-dashed) are displayed as fitted to our fiducial model. Two best fits are applied to the data, 
a linear fit and a lognormal fit, and are shown as purple and red lines. With the exception of the lognormal fit, the slopes above the 
turn-over mass are given in the upper left corner. 
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Fig. 17. Initial mass functions for the models with UV. The UV flux used in these models is 10^^ Gq. Similar to the figures before. 
Figs. [14] to [16] the images display the time averaged IMFs between 1 and 3 free-fall times, where tff = 10^ yr. From left to right, 
the cosmic ray rate increases from 1, 100 to 3000 x Galactic. From top to bottom, the X-ray flux increases from 0, 5.1, 28 to 160 
erg s~' cm^^. In each image, for comparison purposes, the Salpeter IMF (green dashed) and the Chabrier IMF (blue dot-dashed) 
are displayed as fitted to our fiducial model. Two best fits are applied to the data, a linear fit and a lognormal fit, and are shown as 
purple and red lines. With the exception of the lognormal fit, the slopes above the turn-over mass are given in the upper left corner. 
These UV models are simulated for a a black hole mass of Mbh - 10^ Mq. 



